Transport rates and momentum isotropization of gluon matter in 
ultrarelativistic heavy-ion collisions 

Zhe Xu 1 * and Carsten Greiner 1 
1 Institut fur Theoretische Physik, Johann Wolfgang Goethe- Universitdt Frankfurt, 
Max-von-Laue-Str.l, D-60438 Frankfurt am Main, Germany 
[< ; (Dated: March 2007) 

o , 

^ ■ Abstract 

To describe momentum isotropization of gluon matter produced in ultrarelativistic heavy-ion 
\ collisions, the transport rate of gluon drift and the transport collision rates of elastic (gg «-> 



gg) as well as inelastic (gg *-* ggg) perturbative quantum chromodynamics- (pQCD) scattering 



. processes are introduced and calculated within the kinetic parton cascade Boltzmann approach of 

> , 

£C) , multiparton scatterings (BAMPS), which simulates the space-time evolution of partons. We define 

m 

CN ' isotropization as the development of an anisotropic system as it reaches isotropy. The inverse of the 

introduced total transport rate gives the correct time scale of the momentum isotropization. The 
contributions of the various scattering processes to the momentum isotropization can be separated 
into the transport collision rates. In contrast to the transport cross section, the transport collision 



rate has an indirect but correctly implemented relationship with the collision-angle distribution. 
Based on the calculated transport collision rates from BAMPS for central Au+Au collisions at 



Relativistic Heavy Ion Collider energies, we show that pQCD gg *-* ggg bremsstrahlung processes 
isotropize the momentum five times more efficiently than elastic scatterings. The large efficiency 
of the bremsstrahlung stems mainly from its large momentum deflection. Due to kinematics, 
2 — > N (N > 2) production processes allow more particles to become isotropic in momentum space 
and thus kinetically equilibrate more quickly than their back reactions or elastic scatterings. We 
also show that the relaxation time in the relaxation time approximation, which is often used, is 
strongly momentum dependent and thus cannot serve as a global quantity that describes kinetic 
equilibration. 
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I. INTRODUCTION 



It is speculated that the quark gluon plasma (QGP) created in Au+Au collisions at the 
Relativistic Heavy Ion Collider (RHIC) is a strongly coupled liquid [1]. Because of strong 
coupling or rather strong interactions, the QGP fluid has a very small viscosity. However, 
questions regarding the source of the strong coupling and its needed strength to generate 
a quasi-ideal fluid remain unanswered. The necessary condition for the onset of perfect 
hydrodynamical expansion is the achievement of local kinetic equilibrium. Although the 
quarks and gluons produced at RHIC are far from thermal equilibrium, kinetic equilibration 
should occur on a short time scale so that the elliptic flow, v 2 , increases substantially [^[[i,!^]. 
In this article we assume that the strong coupling and thermalization are a consequence of 
frequent collisions among gluons on a semi-classical level. We recently developed a new 
on-shell parton cascade code, BAMPS (Boltzmann approach of multiparton scatterings) Q], 
which is a microscopical relativistic transport model that solves the Boltzmann equation for 
partons that are produced in ultrarelativistic heavy-ion collisions. The included interactions 
can be elastic gg <-> gg processes or inelastic gg <-> ggg pQCD bremsstrahlung processes. 
Although the total perturbative quantum chromodynamics- (pQCD) scattering cross section 

,. — „_,^ a - - 

to generate sufficiently large elliptic flow v 2 [6J. Our goal is to understand the theoretical 
mechanism for the fast equilibration of gluons, which are observed numerically. 

In kinetic theory there are two competing processes that affect kinetic equilibration. The 
first is when particles stream freely between two subsequent collisions. In an expanding sys- 
tem free streaming drives the system out of equilibrium. This is the case in a one-dimensional 
Bjorken expansion, which most likely occurs early on in ultrarelativistic heavy-ion collisions. 
The second one involves collisions that make the particle momentum kinetically isotropic 
and thermal. Here one has to take into account the distribution of collision angle be- 
cause large-angle collisions should contribute more to momentum isotropization. We define 
isotropization as the development of an anisotropic system as it reaches isotropy. A transport 
cross section [?], H] was introduced, either in the form 

a*= [ dO^- sin 2 6 (1) 
do 
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or 

where 9 denotes the collision angle as a pertinent quantity that measures the contributions of 
various collision processes to kinetic equilibration. Although kinetic equilibration is observed 
locally in the comoving frame of the expanding system, the transport cross section is usually 
calculated in the center-of-mass (cm.) frame of individual colliding particles. The changes 
in momenta after the collision appears different in each respective frame. Therefore, the 
transport cross section may not be fully appropriate for characterizing kinetic equilibration. 

A widely used, yet simpler, method to characterize kinetic equilibration is to calculate 
or estimate the relaxation time r re i 0, [h], 11, 12, 13, Q, 15]. In the relaxation time 



approximation the collision term is expressed by (/ cq — /)/T re i, where r re i is assumed to 
be momentum independent and is then a global quantity that characterizes the kinetic 
equilibration time scale. However, the validity of the approximation must be verified. 

In this article we derive a mathematical method of quantifying the contributions of vari- 
ous processes to the momentum isotropization. For this we define the transport rate, which 
is the momentum average of the particle density f(x,p). The particle density is found 
within the parton cascade as a solution of the Boltzmann equation. Moreover, it will be 
shown that the inverse of the total transport rate gives the global time scale of momentum 
isotropization. In Sec. [Ill we mention the operation of the employed parton cascade BAMPS 
and improvements made in it. The initial condition of gluons, as an input for the parton 
cascade, is discussed in Sec. IIHI We show results on thermal equilibration and momentum 
isotropization of gluons in Sec. IIVI for a central Au+Au collision at RHIC (^/i = 200 GeV). 
The inclusion of quarks into the parton cascade is straightforward and the results will be 
presented in another article. In Sec. |V]we define the transport rates, which determine con- 
tributions of various processes to the momentum isotropization and derive their relations to 
the transport cross sections. We present in Sec. IVII the numerical results on the transport 
rates. The transport rate of gluon drift is computed and compared with the one when ob- 
tained assuming Bjorken boost invariance. To show the importance of the bremsstrahlung 
processes in thermal equilibration, we carry out simulations with and without these pro- 
cesses for comparison. The quantitative difference in the momentum isotropization for both 
simulations is manifested by the ratio of the total transport rates. The ratio, which turns 
out to be approximately 6, is used to perform a third type of simulation in which only elastic 



scatterings with artificially enlarged cross sections are included. Although such large cross 
sections are not physical, they verify our main finding: the total transport collision rate is 
the key quantity determining momentum isotropization. Despite the process type as long 
as the total transport collision rate is the same, the momentum isotropization is also the 
same. At the end of Sec. I VI I we demonstrate that the relaxation time approximation is not 
suitable for the quantification of the time scale for kinetic equilibration. A summary of our 
findings is given in Sec. IVHI Detailed expressions for calculating the transport rates are 
derived in Appendix [A] 



II. BAMPS AND SETUP 



The structure of the parton cascade BAMPS is based on the stochastic interpretation 
of the transition rate 0, Q, H, Q- This interpretation ensures that detailed balance is 
not violated, which is nontrivial when the geometrical concept of cross section is used 19] , 
especially for multiple scatterings like ggg «-> gg. BAMPS subdivides space into small cell 
units. In each of which we separately evaluate the transition probabilities of all possible 
gluon pairs and triplets to see if a particular scattering (or transition) occurs. The smaller 
the cells the more local transitions can be realized. However, the smaller cells contain fewer 
particles and thus have larger statistical fluctuations in their calculated transition rates. To 
achieve a high-enough number of pairs and triplets of gluons in a cell, we adopt a test particle 
technique, which amplifies the (pseudo)gluon density by a factor of N test . Accordingly, the 
cross sections have to be reduced by the same factor to obtain the same physical mean free 
path [5]. In this article the transverse length of a cell is a constant of Ax = Ay = 0.25 fm 
and the longitudinal length Az is half of that in 5j, so for a cell of the center of the collision 
Az pa O.lt, where t is the running time of the evolution of gluon matter. N test is set to 280, 
which ensures that there are on average 15 test particles per cell. 

The differential cross section for the elastic pQCD scatterings of gluons is given by 

da" 9 _ 9na* 



Three-body gluonic interactions are described by the effective matrix element 
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| - M — 1 = T (qi + mlY kiKkl - q± y + ml] " ^ (4) 
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where g 2 = 47ra s . a s is set to 0.3 in contrast to the running coupling used in Ref. [5|. and 
h± denote the perpendicular component of the momentum transfer and of the radiated gluon 
momentum in the center-of-mass frame of the collision, respectively, y is the momentum 
rapidity of the radiated gluon in the center-of-mass frame, and A 9 is the mean free path of 
a gluon. 

We regularize the infrared divergences by introducing the Debye screening mass 

m 2 D = 16na s N c J -^~^f 9 (5) 

(N c = 3), which is calculated locally using the current gluon density obtained from BAMPS. 
n g eneral, the Debye screening mass should depend on the direction of the gluon propagator 



221 ] . If the gluon distribution f g significantly deviates from its isotropic shape, the Debye 
screening mass may even become neg ative, which leads to instabilities in certain modes of 



the soft gauge field [23j, |24], |25|, |26|, |27|, |28( . These instabilities and their proper inclusion are 



beyond the scope of the present article. We have simplified the problem by removing the 
directional dependence of the Debye screening mass. 

The suppression of the radiation of soft gluons due to the Landau-Pomeranchuk-Migdal 
(LPM) effect , |2l| is included using the step function in Eq. There the time of the 



emission, ~ -^-coshy, should be smaller than the time interval between two scatterings or 
equivalently the gluon mean free path A g . This leads to a lower cutoff for k± and a decrease 
in the total cross section or the transition probability. 

Compared to the default setup in Ref. j^], further improvements have been made. To 
calculate the Debye screening mass in a local region more accurately, we make use of 
the polar symmetry in central collisions and divide the transverse plane in each Az-bin into 
rings: the first ring has a radial size of < Xt < 1-5 fm (iy being the transverse radius), 
and the following rings have transverse radial widths of 1 fm. The rings are regarded as 
local regions in which the Debye screening mass is evaluated. 

The local collision rates of all interaction channels, the sum of which is the inverse of the 
m ea„ free path that .node.s the LPM effect, we r e eva.uated in Ref. fl in individual ce.,. 
This leads to large fluctuations in the mean free path in cells with few (test) particles. To 
reduce these fluctuations we take the averaged value of the collision rates over all the cells 
within individual rings. In addition, transverse velocities of rings are taken into account to 
calculate the collision rates in the comoving frames. 
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Moreover, we assume that if the energy density, which is calculated locally in the comoving 
frame, sinks below 1 GeV/fm 3 , particles in that region no longer interact, so they propagate 
freely. At this stage a hadronization procedure should be applied, which is planned as a 
future project. 

We concentrate on the central region of the full reaction, which is defined as a cylinder 
with < xt < 1.5 fm and —0.2 < rj < 0.2 where rj denotes the space-time rapidity 
r) = The longitudinal extension of the cylinder is, thus, Az — 2t tanh(0.2) « 0.4 1. 

The parameters for bounding the cylinder are found by balancing between having a small, 
local region and avoiding high statistical fluctuations. Results in this region are obtained 
by averaging over the various ensembles. 



III. INITIAL CONDITIONS 



Initial gluons are taken as an ensemble of minijets with transverse momentum greater 
than 1.4 GeV, which are produced via semihard nucleon-nucleon collisions (29]. Details of 
the distribution of the initial gluons in space and time can be found in Ref. 5|. Using 
Glauber-geometry and assuming independent binary nucleon-nucleon collisions, the gluon 
number is initially about 700 per momentum rapidity. These gluons take about 60% of the 
total given energy entered in a central Au+Au collision. The lower momentum cutoff is taken 
as a parameter to fit the experimentally measured final transverse energy at midrapidity (see 
Fig. ED. 

For simplicity's sake one may assume that the two gold nuclei are extremely Lorentz 
contracted with zero width. Assuming that on-shell gluons are immediately formed (i.e., 
without any formation time) at the same time when the corresponding nucleon-nucleon 
collision occurs, all initial gluons are positioned at z = fm at t = fm/c. Subsequent free 
streaming would immediately order the gluons with the momentum rapidity y to a spatial 
slice with the space-time rapidity rj being equal to y. In the comoving frame of each spatial 
slice gluon momentum has only a transverse component and it has a highly anisotropic 
distribution. 

At RHIC energy each of the colliding gold nuclei has a small but nonvanishing longitudinal 
extension of about 0.2 fm. Therefore, gluons are primarily produced at z = fm (or rj = 0) 
at t ~ 0.1 fm/c when the two nuclei overlap fully. Note that t = fm/c is when two nuclei 
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are just touching. In contrast to the simplified case mentioned above in reality there is a 
significant smearing in the gluonic rj — y correlation for times t < 0.2 fm/c. For instance, 
gluons with y ^ will also appear in the central slice with rj = for a while. The rate of 
smearing disappearance is shown in Fig. [U where the spectra of transverse and longitudinal 
gluon momenta are depicted during initial free streaming in absence of secondary collisions. 




012345678 



FIG. 1: Transverse and longitudinal gluonic momenta spectra during initial free streaming. The 
thick histogram is the distribution of \p x \ at a time of 0.1 fm/c, whereas the thin histograms are 
the distributions of \p z \ at times of 0.1, 0.2, 0.22, and 0.24 fm/c, respectively, from top to bottom. 
The results are obtained from the central region. 

The spectra are obtained in the central region (0 < xt < 1-5 fm and —0.2 < rj < 0.2). We 
see that the \p z \ spectrum changes quite drastically after 0.2 fm/c, the same point when two 
gold nuclei cease to overlap and the production of minijets is completed. The free streaming 
of high \p z \ gluons away from the central region leads to strong, continuous suppression in 
the \p z \ spectrum. Corresponding to this suppression, the changes in the \p x \ spectrum at 
low transverse momentum are, however, tiny compared with their absolute values at low \p x \. 
At high \p x \ the change in time is negligable because there is no initial transverse expansion 
for large nuclei. Therefore, in Fig. [T]the \p x \ spectrum is depicted only at t — 0.1 fm/c. We 
note that at large \p z \ the suppression stops completely when all the particles with high p z 
(or with high y > rjj, = 0.2) have left the small, but finite, rapidity window [—tji, : r)b]- Then, 
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only particles with lower rapidity y remain in the central region. The time a gluon needs to 
leave the central region is, for instance, At = to tanh^/^tanhy — tanh^f,) when the gluon is 
produced at Zq = fm at to- We see that the larger the momentum rapidity y, the smaller 
At. For a gluon with y = 1 and to = 0.1 fm/c it takes At = 0.035 fm/c to leave the central 
region. 

Comparing the particle spectrum of \p z \ with that of \p x \, the momentum distribution is, 
strictly speaking, at no time isotropic during the initial free streaming. The characteristic 
hump in the \p x \ spectrum at 1.4 GeV arises from the requirement that the transverse 

minijets should be greater than 1.4 GeV. Choosing other initial 



momentum of the original 
conditions like in HIJING 
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311 ] or the color glass condensate [32| would change the shape of 



the initial momentum distribution. However, even though the momentum distribution might 
be isotropic during the continuing suppression of high \p z \ gluons, the further suppression 
leads to a deviation in the momentum distribution away from isotropy within a very short 
time of ~ 0.1 fm/c. Therefore, we can conclude that free streaming leads to rj « y regardless 
of the initial rj — y correlation. The gluon momentum distribution after short-time free 
streaming is, in general, neither thermal nor isotropic. 

In this article we introduce an additional formation time [5J for every minijet, Atj = 
coshyArf « coshy • 1/pr, which models the prior off-shell propagation of the gluons to 
be freed in individual nucleon-nucleon collisions, where coshy denotes the Lorentz factor. 
Within At/ we assume that the virtual gluon does not interact and, therefore, moves freely 
at the speed of light. Gluons with large \p z \ have in turn a large Lorentz factor and, thus, 
a large formation time. Although most of these gluons are produced in the central region, 
they are far from the central region when they materialize as on-shell partons because of 
the assumed off-shell propagation. Because we count particles only if they are on-shell, i.e., 
interactive, the initial gluon momentum distribution at 0.2 fm/c differs from that shown 
in Fig. [TJ However, it is practically identical to that at 0.24 fm/c, which is not isotropic. 
When the initial conditions are chosen accordingly and a simulation including the pQCD 
bremsstrahlung processes is performed, we obtain dEx/dy ~ 640 GeV at midrapidity with 
a final time of 5 fm/c, at which the energy density of gluons decreases to the critical value 
of 1 GeV/fm 3 . Our dE^jdy at y = from the simulation is comparable with that found in 



experimental measurements at RHIC 



33| (see Fig. E>- 
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IV. MOMENTUM ISOTROPIZATION AND KINETIC EQUILIBRATION 



Kinetic equilibration is a process in which the particle momentum becomes isotropic and 
thermal, which has an exponential distribution. Momentum isotropization is part of kinetic 
equilibration and is reached before full kinetic equilibrium 
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271 ]. (Strictly speaking, 



full kinetic equilibrium can be achieved only for a static, nonexpanding system.) In this 
article we concentrate on the contribution of collision processes to momentum isotropization 
and kinetic equilibration of gluon matter in ultrarelativistic heavy-ion collisions. 

As demonstrated in Sec. IHIt the initial free streaming (or the off-shell propagation) of 
gluons with high momentum rapidity y makes the momentum distribution anisotropic, even 
if it appears momentarily isotropic. Initially in the central region most gluons move in the 
transverse direction. Secondary collision processes gradually force them into the longitu- 
dinal direction, which gives a positive contribution to momentum isotropization. However, 
whenever a gluon switches to the longitudinal direction, its momentum rapidity grows and 
the gluon tends to drift out of the central region. This gives a negative contribution to mo- 
mentum isotropization in the local region. Although gluons with the same (regardless of ± 
sign) momentum rapidity drift from their neighboring slices into the central slice, this cannot 
completely compensate for the loss in the central region. The reason is that thermalization 
occurs earlier in the central slice than in the outwards regions corresponding to Bjorken's 
picture of boost-invariance in the space-time evolution of the parton system [341]. In the 
transverse direction, however, there is no transverse flow at the beginning of the expansion. 
Therefore, no net drift of gluons occurs in the transverse direction. The difference in the 
gluon drift in the longitudinal and transverse directions leads to a situation in which the net 
effect of the drift has a negative contribution to momentum isotropization and the stronger 
the momentum isotropization, the larger the negative contribution of the particle drift. At 
later times, when three-dimensional expansion takes place, there is also a net particle drift 
in the transverse direction and the negative contribution of the particle drift to momentum 
isotropization decreases. 

In this section we first demonstrate momentum isotropization and kinetic equilibration of 
gluons in central Au+Au collisions at RHIC energy. The various contributions of collisions 
and drift to momentum isotropization will be analysed in detail in the next section. Figure 
[2] depicts the transverse and longitudinal gluon momenta distributions in the central region 
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at four different times throughout the evolution of gluon matter. The results are obtained 
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FIG. 2: Transverse (thick histograms) and longitudinal (thin histograms) gluon momentum spectra 
in the central region at various times throughout the evolution of gluon matter. Results are obtained 
from the simulation when both elastic and inelastic pQCD-scattering processes are included. 



when both elastic and inelastic pQCD-based scattering processes are included. We first see 
that the momentum distribution continuously isotropizes and thermalizes over time. Due 
to the expansion full thermal equilibrium cannot be achieved if the collision rate is finite. A 
certain mismatch in p z and p x must exist due to the counteraction between the expansion 
and the collisions. 

Second, Fig. [2] shows that there is an exponential distribution before the system becomes 
isotropic. It is almost impossible to distinguish the gluon momentum isotropy time scale 
from the thermal time scale. It seems that when collisions drive the particle momentum 
close to isotropy, the momentum distribution is already practically thermalized. In general, 
momentum isotropization happens on a shorter time scale than kinetic equilibration. The 
difference in the time scales of both dynamical processes depends on the initial condition 
for gluons. 

Kinetic equilibration for the softer gluons is completed earlier than that for the harder 
gluons. This is obvious for elastic gg <-» gg scattering processes because the momentum 
transfer in collisions is typically the Debye screening mass. Therefore, the hard gluons 
cannot be deflected as strongly as the soft gluons. However, in the inelastic pQCD gg <-> 
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ggg collisions, which we will prove are the dominant processes in kinetic equilibration, the 
difference in the momentum degradation for soft and hard gluons is miniscule due to the 
production or absorbtion of an additional gluon. Averaging the various kinetic equilibrium 
times for soft and hard gluons the momentum distribution becomes isotropic and thermal 
at 1 — 2 fm/c. Furthermore, we clearly see that the distributions become steeper with time, 
which indicates the ongoing cooling of the system related to quasihydrodynamical behavior 
due to the subsequent work done by the expanding system. 

To understand the role of the inelastic pQCD gg <-> ggg processes in kinetic equilibration, 
we also carry out calculations in which gluons interact only via elastic scatterings. The initial 
conditions are the same as those when inelastic collisions are included. The results are shown 
in Fig. [31 which has the same structure as Fig. [2J The difference in the results depicted 




Px/z [GeV] P A [GeV] 

FIG. 3: Transverse (thick histograms) and longitudinal (thin histograms) gluon momentum spectra 
in the central region at various times throughout the evolution of gluon matter. Results are obtained 
using elastic only pQCD-scattering processes. 

in both figures can be immediately seen. The spectra in Fig. [3] only show a small change 
throughout the entire evolution of the system and are still highly anisotropic and are not 
thermalized as late as 4 fm/c. The evolution resembles that of free streaming. 

The kinetic equilibration time strongly depends on whether the pQCD bremsstrahlung 
processes and their back reactions are taken into account. The pQCD bremsstrahlung 
processes and their back reactions isotropize the momentum more efficiently than elastic 
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collisions and, thus, play an essential role in early thermalization of gluons in heavy-ion 
collisions at RHIC. As seen in Fig. [U the pQCD cross section of gg — > ggg processes, 
including LPM suppression (dashed curve), is smaller than that of elastic scatterings (solid 
curve) and much smaller than the cross section obtained in the simulation with elastic-only 
scattering processes (dotted curve). Although particle production in inelastic processes can 
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FIG. 4: Time evolution of pQCD cross sections. The solid and dashed curve show the pQCD 
cross section for gg — > gg and gg — > ggg collisions, respectively, when both elastic and inelastic 
scattering processes are considered. The dotted curve gives the cross section of gg — > gg colli- 
sions in the simulation, including elastic-only pQCD scattering processes, whereas the dash-dotted 
curve indicates the cross section (divided by a factor of 6) of gg —* gg collisions when elastic-only 
scattering processes with artificially large cross sections are included. 



enhance the number of collision centers and, thus, effectively shorten the mean free path of 
particles, chemical equilibration will balance the production by the annihilation of particles 
to avoid oversaturation. The fact that the cross section of the pQCD bremsstrahlung process 
is small, but its kinetic equilibration efficiency is large, demonstrates that cross sections or 
collision rates are not the correct quantities to describe the contributions of various processes 
to kinetic equilibration. The collision- angle distribution must be at least taken into account. 
Defining the correct quantity is one of the main purposes for this article. 

In Fig. H] the large difference in the total cross sections of elastic scatterings for various 
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simulations is shown. Because 

(vgg-vs) ~ ml(l + 47711/3)' (6) 

the difference in the cross sections arises from the difference in the development of the Debye 
screening mass in the various simulations, is calculated dynamically according to 



(J5]) and, thus, is roughly proportional to Jn/(p) = n/y/e, where n and e are the number 
and energy density of gluons, respectively. 

We consider two extreme cases of expansion with initial conditions that possess the lon- 
gitudinal boost invariance. One case is free streaming, for which n as well as e decrease as 
r _1 , where r = \Jt 2 — z 2 is the proper time. Thus, ttid decreases as r -1 / 2 . In the other 
case of a one-dimensional ideal hydrodynamical expansion, n decreases as r -1 , whereas e 
decreases as r~ 4//3 . Therefore, tud decreases as r" 1 / 3 . In a viscous hydrodynamical ex- 
pansion the decrease of over time falls between the two cases. The time dependence 
of the Debye screening mass in a real expansion starting out of thermal equilibrium and 
undergoing thermalization is more complicated. Whereas kinetic equilibration drives the 
density distribution of gluons to its thermalized shape, which affects the calculation of 
(JSJ), chemical equilibration, which is not taken into account above, will enhance or reduce 
the gluon number, which in turn enhances or reduces mp. 

Figure [5] shows the time evolution of the Debye screening mass in various simulations. 
The results are obtained in the central region where t ~ r. All the curves in Fig. [5] decrease 



with time. Similar calculations for the Debye screening mass have been done in Refs. [351.136] 
employing the parton cascade VNI/BMS. From Fig. [5]we see that the evolution depends on 
the simulation type: in the simulation including both elastic and inelastic pQCD scatterings 
(solid curve) mo decreases slower than t -1 / 3 due to gluon production in the course of chemical 
equilibration; in the simulation with elastic-only pQCD collisions the decrease of mo (dotted 
curve) is slightly stronger than t -1 / 2 , which indicates again that the expansion of gluons in 
this simulation resembles that of free streaming; the third simulation includes elastic-only 
scatterings with artificially large cross sections and shows the same kinetic equilibration as 
that in the simulation including both elastic and inelastic pQCD collisions (see Fig. [6]). The 
Debye screening mass in this simulation (dash-dotted curve) decreases between t~ 1 / 2 and 

Returning to the kinetic equilibration analysis, the gluon kinetic equilibration time can 
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FIG. 5: Time evolution of the Debye screening mass. Results are obtained from the simulation 
with both pQCD elastic and inelastic collisions (solid curve), with elastic-only pQCD collisions 
(dotted curve), and with elastic-only collisions using large cross sections (dash-dotted curve). 

in principle be determined quantitatively by studing the entropy production. Because the 
entropy can be hardly extracted from any microscopic cascade, we concentrate on momentum 
isotropization of gluons. Choosing minijets production as the initial condition, momentum 
isotropization and kinetic equilibration time scales are almost identical (as seen in Fig. [2]). 

To quantify momentum isotropization we have to choose an appropriate momentum- 
distribution moment Q. For instance, Q := (p 2 z / E 2 ) is used to describe momentum 
isotropization. Later we briefly discuss the consequences of Q = (\p z \/E), to see how 
sensitive the results are to different descriptions of momentum isotropization. In Fig. [6] mo- 
mentum isotropization with Q = (p 2 /E 2 ) is depicted. The average is taken over all gluons 
in the central region. In Fig. [6] we see that Q relaxes its equilibrium value of 1/3 when the 
inelastic processes are included, whereas it still deviates from its equilibrium value at the 
time 4.5 fm/c when only elastic pQCD scatterings are considered. These results agree with 
the momentum spectra time evolution shown in Figs. [2] and [3j The dash-dotted curve in 
Fig. [6] depicts the momentum isotropization considering elastic-only collisions with artifi- 
cially large cross sections and is almost the same as the solid curve, which implies they have 
basically the same kinetic equilibration. The third simulation is detailed in Sec. IVII 
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FIG. 6: Momentum isotropization. Results are obtained from the simulation with both elastic 
and inelastic pQCD-scattering processes (solid curve), with elastic-only pQCD-scattering processes 
(dotted curve) and with elastic-only scattering processes using artificially large cross sections (dash- 
dotted curve). 



The momentum isotropization fit is found using the relaxation formula 



F<*) = 5 + 



exp 



t-t 

#rel(£() 



(7) 



F(t) is equal to Q(t) only at t = to- For every fixed to the relaxation time 8 re \ is constant 
with respect to t. Using 9 re \ = 0.9 fm/c at to = 0.3 fm/c up to 1.0 fm/c and 9 re \ = 2.4 fm/c 
at to = 1.2 fm/c for the rest, F(t) is a perfect fit for the solid curve in Fig. [6] An isotropy is 
achieved at about 1.0 fm/c in the simulation that includes both elastic and inelastic pQCD- 
scattering processes. This time scale is consistent with that extracted from the momentum 
distribution (see Fig. [2]). Within our parton cascade description early thermalization occurs 
at roughly 1 fm/c for the initially nonequilibrated gluon matter at RHIC. 

The relaxation time # re i is generally time dependent. Because a local fit requires that the 
time derivatives of F(t) and Q(t) are equal at t = to, which leads to 

= H/) = -KJI/n) - (),,,) — J 

t=to 



Q(t) 



t=to v-v-uy ^e iel (t 

where Q eq = 1/3, 6 re \ can be calculated as [changing t to t in Eq. 

<m _ i 

Qeq-Q(t) e TcX {t)' 



(8) 



(9) 
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Equation (Q expresses the relaxation rate of momentum isotropization l/6 re \ as a function 
of time, which will be separated analytically into terms corresponding to the particle drift 
and the various scattering processes. We will derive the so-called transport rates, which 
precisely quantify the contributions of various processes to momentum isotropization. 



V. TRANSPORT RATE 

To introduce Q at a certain space point £ one has to consider its comoving frame. For 
the coordinate x and the momentum four-vector p in the £'s comoving frame, Q is defined 
by 



Q(t) := {§,] 
where the local number density is 



1 1 >r ' p ni mt, P )\ s=0 , (io) 



2 = 



nJ (2vr) 3 £ 2 



«(*)=/ ^3 f&*>P)\**' (11) 

In practice Q is evaluated within a volume element, which is small compared to the volume of 
the expanding system but is large enough so that it still contains a large number of particles. 
For the calculations shown in Fig. [6] we used the central region bounded by xj< < r& = 1.5 
fm and \rj\ < rjb = 0.2. Correspondingly ( Tl0|) and ( ITTi) must also be adjusted and detailed 
expressions are derived explicitly in Appendix [A] As a simplification we now consider the 
limit r b — > and r\ h — > where definitions ffTU]) and ffTT]) can be used. 
Taking the time derivative of Q(t) yields 



1 y d 3 p df 



(12) 



x=0 



We replace df/dt in Eq. ffT2]) by 

^ = -| • V/ + C 22 + C 23 + C 32 (13) 

from the Boltzmann equation, where — % ■ V/ corresponds to particle drift and C 22 , C 2 3, 
and C 32 denote the collision terms corresponding to gg — > gg, gg —> ggg, and ggg — > g>^, 
respectively. It is obvious that the contribution of the various processes to Q(t) is additive. 
We rewrite Eq. ([12]) 

Q(t) = W dTi{t (t) + W 22 (t) + W 23 (t) + W w (t) , (14) 
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where PVdrift, W22, W23, and W& correspond to particle drift, gg — > gg, (7(7 — > ggg, and 
fl'fl'fl 1 "~ > collision processes, respectively. According to Eq. (jUJ) we obtain 

1 



where we define 



RlAt) + R'm + i£(f) + i£(f) , (15) 



for i = drift, 22, 23, and 32. One sees that the relaxation rate of momentum isotropization 
1/6 TC \ is separated into additive parts corresponding to the particle drift and the various 
collision processes. -Rdrift i s called the transport rate of particle drift, whereas R 22 , -R23, and 
-R32 stand for the transport collision rates of their respective interactions. Extending this 
to more than three-body processes is straightforward because the collision term is additive. 
We note that R\ r (shown below) depends on the definition of Q. When one changes Q from 
Q = (pI/E 2 ) to Q = (\p z \/E), the form of Rf changes accordingly. 



A. R* v 



drift 



Except for static systems the drift term in the Boltzmann equation fll3p generally con- 
tributes to Q{t). Wdrift is given by 

^2 



^ftW^/^l-v; 



(17) 



n J (2tt) 3 E 

Assuming Bjorken's space-time picture of a central ultrarelativistic heavy-ion collision |34|, 
we can use the relation 

l.Vf* P - d -l = - P -^- (18) 
E E dz t dp z y ' 

found in Ref. [9j. Inserting Eq. (TI51) into ( 1171) and performing partial integrals we obtain 



-^drift(*) 



Q(t) 



E A 



(t) 



(19) 



[Qeq - Q(t)] t _ 

Equation (jT9l) shows that -Rjfrift i s negative, which agrees with our conclusion in the previous 
section. Using the approximation (p A /E A ) « Q 2 , we see that the larger the Q, the larger the 



ptr 
"-"-drift- 
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B. R tr 



22 



Changing p to p±, W22 becomes 

w -® = kf$? 1 k c »- (20) 

where C22 does not contribute to the second integral in Eq. (fl2l due to particle number 
conservation in elastic collisions. The same holds for the sum of C23 and C32 in chemical 
equilibrium. Inserting the explicit expression of the collision term 

° 22 = Wj ^4 / «/i/2l-^i'2'^i2l 2 (2 7 r) 4 5 (4) (p'i +P2 -Pi -ft) 

"2X / dT 2fif2^JdT[dT' 2 \M 12 ^ 2/ \\27r)^\p 1+ p2-p[-p^ (21) 

(dTi = d 3 pi/(27r) 3 2Ei for short) into Eq. gives two terms, which indicate the "gain" 
and "loss" in momentum isotropization. 
The loss term is 

n2 



^jdT 1 dT 2 f 1 f 2 ^2sa 22 = nU e ^a 22 \ , (22) 



where 



^22 := ~ / « |.M 12 ^ 2 f (2tt)V 4 )( Pi +J9 2 - p[ -p' 2 ) (23) 

is the total cross section, s is the invariant mass of the colliding system, v re \ = s/2EiE 2 is 
the relative velocity, and ( ) 2 symbolizes an ensemble average over incoming particle pairs. 
In BAMPS f(x,p) = Ysi^lx ~ Xi(t)}6^ (p — pi) and we evaluate the averages ( ) 2 in local 
cells by running over all particle pairs in the cells. Each cell has a small volume to ensure 
local collisions and has a sufficient number of (test) particles to achieve adequate statistics. 
The iy 22 's gain term is n(v Te \a22)2, where 

Z22 ■= ~ J « |f |AW 2 f (2vr)V 4 )( Pl + P 2 -p[ -p 2 ) , (24) 

which, like Eq. ( 123]) . is an integral over all possible states of outgoing particles. Equation 
( |24|) was obtained by exchanging the primed and unprimed variables in Eq. fT2T|) . Except 
for p'^ z /E'i all variables and functions in Eq. ff2~l"j) are Lorentz invariant. Particularly we 
find dV' x = dY'{ = d 3 p{ / (27r) 3 2£i* = dQ*dE[* E'{ /2(2tt) 3 , where p[* is the four-momentum 
of an outgoing particle in the center-of-mass frame manifested by pi and p 2 of the incoming 
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particles, and Q* denotes the solid angle relative to the collision axis in the center-of-mass 
frame. Integrating over dY' 2 = dT' 2 * using the four- dimensional 5 function gives 



*»= /dn'^^f, (25) 



■. df22 Pu 

dtt* E'i 1 

where p lz and E[ are the Lorentz transformed quantities from p* and, thus, functions of 
s, and p. The (3 = {p\ +P2)/(Ei + E 2 ) denotes the relative velocity of the center-of-mass 
frame of colliding particles to the laboratory frame where Q is defined. 
We finally obtain 

Q^ro^m ("(-/^|> 2 - ! ("4->J ■ <*> 

where the momentum isotropization gain and loss terms are clearly seen. The relationship 
to the collision- angle distribution is implicitly contained in Rjg 2 - When the collision rate is 
defined as 

R 22 = n{v tc ia22)2 , (27) 

then we call R^ 2 the transport collision rate of elastic scatterings. 

The transport collision rate R 22 in Eq. fl26|) . in general, differs from n(v re \0~22)2, where 
(T22 is defined in Eq. (JTT ) or They match only if the laboratory frame is identical to 
the center-of-mass frame of colliding particles. To demonstrate this we consider the special 
case in which half of the particles move along the positive z axis and the other half of the 
particles move along the negative z axis and all the particles have the same energy E such 
that 

f(x,p) ex 5(p x )5(p y )5(p z - E) + 5( Px )5(Py)S(Pz + E) . (28) 

In this case the laboratory frame is the same as the center-of-mass frame for every colliding 
pair. Thus p^jE'l = cos 2 8* and p\ z jE\ = 1. We then have 

i4 = ^(iw4> 2 (29) 

where a 22 is given in Eq. (0Q). It is easy to verify that Eq. ( 1291) does not depend on the 
direction of the initial momentum. The only necessary conditions are that all the particles 
move along the same (regardless of ± sign) direction and have the same energy. Also, if 
Q — (\Pz\/E), R22 will be changed to 

R^ = 2n(v rcl a t 2 I 2 ) 2 , (30) 
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where <y\\ is given in Eq. (j2j). 

The reason R^ 2 is called as the transport collision rate now becomes obvious because Eq. 
( I26D is the generalization of the simplified formula na tr , which is referred to in the literature 

nn 

as the transport collision rate [7, 8]. 

To understand the physical meaning of the transport collision rate it is reasonable to 
interpret as the rate per particle at which particles experience elastic collisions to be- 
come isotropically distributed in momentum space, because R 22 contributes to momentum 
isotropization according to Eq. ([15!) . For ultrarelativistic particles the inverse of R^ is the 
mean path (or time) that particles should travel to become isotropic, and R22/R22 is the 
average number of collisions, which each particle needs to drive the particle system into 
isotropy in momentum space. 

To confirm this interpretation we calculate -R22 assuming that the collision angle is isotrop- 
ically distributed. We then obtain R% 2 = R22 via Eq. (1291) or (1301) for the special case (1281) . 
This indicates that each particle needs only one collision to drive the particle system into 
isotropy in momentum space if the distribution of the collision angle is isotropic. A more 
general case occurs during equilibration. The energy spectrum of particles tends to be a 
Boltzmann distribution. Rarely found high-energy particles need on the average more than 
one collision to become isotropic, even if the distribution of the collision angle is isotropic. 
The reason is that a particle with high energy always collides with low-energy particles. 
The relative velocity of the center-of-mass frame to the laboratory frame is large and, thus, 
the Lorentz boost has a strong effect. In the laboratory frame deflection in the momentum 
of high-energy particles is narrower in the forward direction. However, low-energy parti- 
cles move perpendicularly to their initial direction and, thus, their momentum deflection is 
large. The averaged effect of the Lorentz boost on the momentum isotropization is, however, 
nontrivial and must be calculated numerically. 

The above hinges on the assumption that the system is static. Expanding systems are 
more complicated because particles flow. Collisions not only deflect the particle momenta but 
also force particles to flow. When we include the flow, which is the particle drift contribution 
to momentum isotropization [see Eq. ffl5l) ]. the momentum degradation of flowing particles 
toward isotropy is slower (~ 6 TC {) than the inverse of the transport collision rate, because 
the transport rate of particle drift is negative in an expanding system. 
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C. Eg, and R% 



Compared with W22 in Eq. (!20|) . W23 has an additional term due to particle production 

^^UM c *-<J^ c - (3i) 

Inserting the explicit formula 

c 23 = / «r 3 l J i^ lV ^ 123 i 2 (27r) 4 5( 4 )(p' 1 + P ; - Pl - P2 - P3 ) 

I ^ 2/l/ 4 / <«i^l2-,l'2'3'| 2 (27T)V 4 )( Pl + P2 - P ; - pi -,4) (32) 

into Eq. (13TT) . we obtain 

3 / - \ / pL \ 1 

- n{V Tel 0-23/2 - Tl &23 J ~ 2 

where 



W 23 (t) = - n(v rc i cr 23 ) 2 - n ( ^rci^f o- 23 ) - - Q(t) n{v ie i cr 23 ) 2 , (33) 



11 2 

0^23 := I «rfr; |f |-M 12 ^ 3 f (2tt)V 4 )( Pi + P 2 -ft -ft -ft) . (34) 

The formula for o"2 3 is just Eq. ( I34p . excluding p'^/E'^. The first two terms on the right- 
hand side of Eq. (l33j) . the sum of which is equal to the first term on the right-hand side 
of Eq. fl3T|) . have similar forms as those in Eq. fl26|) [multiplying Q cq — Q(t)} for W 2 2- The 
coefficients for the momentum isotropization gain and loss terms, 3/2 and 1, indicate that 
in a 2 — > 3 collision the ratio of the gained to the lost particle number is 3/2. The last term 
in Eq. (I33I) stems from pure particle production. The coefficient for this term, 1/2, comes 
from the sum of the gain and loss terms in the particle production process. For a general 
M — > N collision the coefficients will be N/M, 1, and (N — M)/M, respectively. Assuming 
that 

(v Iel ^a 2 ^j ~ (^^j (Vre\ ^23} 2 ^ Q (t) (v lel a 23 ) 2 (35) 

and then comparing W23 in Eq. (|33l) to W22 in Eq. ( 1261) [multiplying Q cq — Q(t)] we realize 
that a gg — > ggg collision is a factor of 3/2 more efficient for momentum isotropization than 
a gg — > gg collision, when o" 2 2 = o~2 3 and 0-22 = cr 23 . The physical reason is obvious: a 2 — ► 3 
collision brings one more particle toward isotropy than a 2 — > 2 collision. 

For the special distribution function ( l28l) we find a relation between the transport collision 
rate and the transport cross section flTJ 

#23 = ^^rclof 3 > 2 . (36) 
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For scattering processes with isotropically distributed collision angles one obtains R\ 3 — 
| R 23 where 

-R23 = n(v rc icr 23 ) 2 (37) 

denotes the collision rate for a gluon undergoing gg — > ggg collisions. Bremsstrahlung effec- 
tively shortens the mean transport path of particles that are becoming isotropic in momentum 
space. Generally, in a 2 — > N process 

N 3 

RIn = J ^n(v iel a^ N ) 2 (38) 

and the larger the number N, the stronger the effect. 

The final expression for W 32 (intermediate steps are analogous to those for W 23 , and C 32 
is found in Ref. [5j) is given by 



1 2 / J 32 \ 1 2 I Viz J 32 \ 1 2 / J : 



W 32 (t) = -n 2 (^^ ) -- n 2 (^§^^ ) +lQ(t)n 2 {— ^ ) , (39) 

3 \8EiE 2 Es 1 3 2 \ E\ 8EiE 2 E 3 j g 6 \ 8EiE 2 E 3 / 3 

where 

1 '2 

^32 := 2j / « |§ |-M 123 ^ lV I 2 (2vr) 4 5 (4) ( Pl + p 2 + - A ~ P 2 ) • (40) 

-Z32 is just Eq. (HOI) , excluding p 2 z /E' 2 . ( ) 3 denotes an ensemble average over triplets of 
incoming particles. 

Comparing W 2 s to VF32, we see that the sum of the last terms in Eqs. (133]) and (|39|) 
originates from the second term in Eq. f[T2"j) but substituting C23 + C32 in for df jdt and it 
should be zero at chemical equilibrium. We obtain 

Tl{v iel <T 23 >2 = \n 2 ( ^ ) (41) 

or, equivalently, -R23 = I-R32, where 

^(m),- (42) 

From Eq. f j4Tl) we derived the collision rate of a gluon experiencing ggg — > (7(7 collisions. 
Assuming that 



we finally have 



E\8EiE 2 E 3 ) ^ ^ { 8EiE 2 E 3 ^ 3 ^ ^ 



3 

W 23 (£) ~ -(^(^rcl 5-23)2 - Q(t)^(^rel 0"23)2) 

32 ^ ^ 3 ^ 8EiE 2 E 3 ^ 3 ^ ^3 \ 8EiE 2 E 3 ^ 3 
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The expansion together with Eq. (HTj) leads to W23 ~ ^W 32 and -R23 ~ §-^32 f° r chemical 
equilibrium. Thus, a 2 — > 3 process should contribute more to kinetic equilibration than 
a 3 — > 2 process because it brings one more particle toward isotropy. If the system is 
out of chemical equilibrium, one expects R^ 3 ~ §^^321 wnere the gluon fugacity X g = 
1 at chemical equilibrium. In an undersaturated system (X g < 1), for instance, particle 
production dominates and, therefore, -R23 i s m uch larger than R 32 - 

For the special case (128]) there is a direct relation between the transport collision rate 
and transport cross section [see Eqs. (I29I) and (|36|) ]. The same should be self-evident for 
-R32 when detailed balance is considered 

R%™\\ ^ = \\n(v jel a%) 2 . (44) 

If the distribution of the collision angle is isotropic, R^ 2 ~ A 9 i? 23 = I-R32, where \ g = 
2R 32 /3R 23 is used. For a iV -> 2 collision (JV > 2) 

2 3 

^2 ~ ^ K R2N = 2 A 9 ^(^el4 r iV>2 , (45) 

which is not proportional to iV in contrast to R 2N in Eq. ( 1381) . For large a 2 — > N process 
kinetically equilibrates significantly more efficiently 

We summarize the main findings derived in this section: 

(i) In Eq. ( fl5i) we showed that the relaxation rate of momentum isotropization is a sum 
of the transport rate of particle drift and the transport collision rates of the various 
scattering processes. 

(ii) The transport rate of particle drift is negative for an expanding medium, which means 
that the particle drift counteracts the momentum isotropization. 

(iii) The transport collision rates of the various interactions found in Eqs. (1261) . ( 1331) . and 
( |39l) [over Q eq — Q(t)} have indirect but correctly implemented relationships with the 
collision-angle distributions. 

(iv) 2 — > N (N > 2) processes isotropize the momentum more efficiently than elastic 
collisions or annihilation processes because the production process brings more than 
two particles toward isotropy in momentum space. 
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(v) The relations between the transport collision rates and the transport cross sections for 
the special case in (1281) 

3 3 3 3 

#22 = 2 n (^cl4 2 )2 , #23 n ( V ™l°2%h , #32 « £ ^ ^^S^ (46) 

are found as long as Q = {p 2 /E 2 ) and the transport cross section is defined by Eq. 
(CQ). However, if Q = (\p z \/E) and the transport cross section is defined by Eq. (T5]), 

3 

R% = 2 n(v iel af 2 ) 2 , R$ 3 = - 2n(v iel a%) 2 , R% « 2 A g n(w rel 4 r 3 ) 2 . (47) 

For the isotropic distribution of the collision angle we find 

3 2 

#22 — #22 , #23 — 2 #23 > #32 ~ 3 #32 • (48) 

VI. RESULTS FROM THE PARTON CASCADE CALCULATIONS 

In this section we present results on the gluon transport rates in the central region of the 
expansion simulated by BAMPS. We then compare the transport rates with those obtained 
from the standard concept of the transport cross sections. The potential dependence of the 
relaxation time on momentum is determined. 



A. Transport rate 

In Fig. 0— -Rjf rift is shown. It cannot be computed by Eq. (Tl7|l because of strong numerical 
uncertainties in calculating V/. Instead, it is obtained by summing ±[Q(t) —p 2 z /E 2 ] over the 
particles, which come into (+ sign) as well as leave (- sign) the central region within a time 
interval of 0.1 fm/c. This causes a large statistical fluctuation. Comparing Fig. [6] to Fig. [7] 
we realize that the larger the Q = (pl/E 2 ), the larger is the (negative) effect of the particle 
drift on momentum isotropization. This confirms our qualitative understanding outlined in 
Sec. HYl The dotted curves estimate the transport rate according to Eq. fTl9|) assuming a 
one-dimensional Bjorken boost-invariance expansion. Q(t) and (p A z /E A ){t) come from the 
parton cascade. At intermediate times our estimates nicely match the numerical results, 
which indicates that the expansion follows a Bjorken expansion. Early in the expansion 
the particle drift is stronger due to free streaming caused by the initial conditions. Later 
on the expansion becomes three-dimensional and particles begin to flow outward in the 
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FIG. 7: Particle drift transport rate (multiplied by —1) from both elastic and inelastic pQCD- 
scattering processes (solid curve), from elastic-only pQCD-scattering processes (dashed curve), and 
from estimates in Eq. (fTUj) (dotted curves). 

transverse direction. The transverse drift of particles with large px is then similar to the 
longitudinal drift of particles with large p z ~ px- The net effect of the particle drift on 
momentum isotropization in a three-dimensional expansion diminishes in comparison to 
a purely longitudinal expansion, as demonstrated in Fig. [7] with the comparison of the 
numerical results to the estimations. 

The numerical results for the transport collision rates are calculated using the expressions 
f l26|) . fl33|) . and fl39|) {the last two are divided by [Q eq — Q(t)]}, and are shown in Fig. [HJ One 
realizes the dominance of the inelastic collisions in momentum isotropization by computing 
the ratio (R% 3 + -R32V-R22' which is about 5 throughout the entire evolution of the system. 
The ratio of to R^ 2 is always larger than 3/2 but nears 3/2 late in the expansion. 
According to Eq. (jUJ) the system is undersaturated early on and eventually reaches chemical 
equilibrium. When we compare -R^s obtained from the various simulations, the difference 
is small, unlike for the cross sections shown in Fig. HI The reason lies in the difference in 
the evolution of the Debye sceening mass for the various simulations as shown in Fig. [5j A 
smaller Debye screening mass leads to a larger cross section but also a smaller collision angle. 
The former causes more frequent collisions and, thus, speeds up equilibration, whereas the 
latter causes inefficient momentum deflection and, thus, slows equilibration. Both contribute 
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FIG. 8: Transport collision rates (R22, -^23 an d R32) from both elastic and inelastic pQCD- 
scattering processes (thick solid, thick dashed, and thick dotted curve, respectively), RW from 
elastic-only pQCD scatterings (thin solid curve) and R22/6 from elastic-only processes using arti- 
ficially large cross sections (thin dash-dotted curve). 

to the transport collision rate so that it is not particularly sensitive to the Debye screening 
mass unlike the total cross section. In Fig. H]the total cross sections for elastic collisions differ 
by a factor of 4 — 6 between elastic-only scatterings and those that include bremsstrahlung 
processes, whereas the corresponding transport collision rates in Fig. [8] are nearly identical. 

The ratios of elastic+inelastic scatterings to elastic-only collisions for the total transport 
collision rate and the transport rate of particle drift are almost identical: the ratio increases 
from 4 at 0.3 fm/c to 9 at 4.5 fm/c. The inverse of the ratio of the momentum isotropization 
time scales in the two simulations is also the same (see Fig. [TT]) . 

Because the change in particle drift is a consequence of particle collisions, one may expect 
that the momentum isotropization is dependent only on the total transport collision rate. 
Gluon kinetic equilibration would always look the same, if the total transport collision rate 
in every evolution was the same at every space-time point. The types of collision processes 
are not relevant, although they are interesting in their own right. We have already shown 
two examples of evolution of gluons in a central Au+Au collision at RHIC energy. The total 
transport collision rate becomes on average a factor of 6 larger if pQCD bremsstrahlung 
processes are included. For another evolution to have the same total transport collision rate 
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as that obtained when pQCD bremsstrahlung processes are included, elastic-only scattering 
processes with larger cross sections, namely da^ildt — 6 do^ GVl /di, were used. If the elastic 
pQCD cross sections obtained from the new simulation were the same as those from elastic- 
only scatterings with pQCD cross sections (see the dotted curve in FigH]), 6 would be an 
appropriate prefactor. The dash-dotted curve in Fig. [4] shows the elastic pQCD cross section 
calculated from the new simulation, which is a factor of 2 smaller than the dotted curve. 
Recalling that the cross sections are dependent on the development of the Debye screening 
mass and that the gluon evolution resembles free streaming for elastic-only pQCD scatterings 
when the Debye screening mass decreases as ~ t _1//2 , artificially large cross sections decrease 
the Debye screening mass from t -1 / 2 to t^ 1 ^ 3 (see the dash-dotted curve in Fig. [5]), which 
implies that the evolution of gluons for large cross sections is a hydrodynamical expansion 
with a finite viscosity. 

-R22/6 for large cross sections is depicted in Fig. [8] and is nearly the same as the transport 
collision rate for standard pQCD cross sections, which proves that the transport collision 
rate for elastic pQCD scatterings is not sensitive to the Debye screening mass. Therefore, 
the total transport collision rates for elastic-only collisions with large cross sections and 
for both elastic and inelastic pQCD scatterings are nearly the same, which implies the 
same momentum isotropization in both simulations. Comparing the time evolution of the 
momentum isotropization (the solid versus the dash-dotted curve in Fig. [6]), we realize that 
momentum isotropization is indeed nearly the same. However, the total cross sections are 
very different (see Fig. H|). At 4.5 fm/c, for instance, (1722)2 + (023)2 ~ 4 mb for elastic and 
inelastic scatterings, whereas (022)2 ~ 45 mb for elastic-only collisions, which is a factor of 
12 larger! 

Because kinetic equilibration and pressure buildup are related we also expect that pres- 
sure buildup does not depend on the type of interactions. Figure [9] shows the time evolution 
of the transverse energy per unit momentum rapidity at midrapidity. The decrease in the 
transverse energy indicates that mechanical work has been done by pressure gradients, which 
are built up during kinetic equilibration. From Fig. [9] one realizes that the time evolution of 
dE^jdy\ y= Q obtained from elastic and inelastic scatterings and from elastic-only scatterings 
with large cross sections are almost identical. This indicates that the ongoing kinetic equi- 
libration and the pressure gradients buildup are the same not only at the collision center as 
already shown in Fig. [6] but also at the central slice of the expansion. There only the total 
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FIG. 9: Time evolution of the transverse energy per unit momentum rapidity at midrapidity 
for both elastic and inelastic pQCD scatterings (solid curve), for elastic-only pQCD scatterings 
(dotted curve), for elastic-only scatterings with large cross sections (dash-dotted curve), and the 
ideal hydrolimit, cIEt / dy\ y= o ~ t -1 / 3 (thin solid curve). 

transport collision rate matters, not the detail of the interactions. 

Whereas the decrease in dEj</dy\ y= o for elastic-only pQCD collisions is very weak, which 
implies slow momentum isotropization, the decrease in dET/dy\ y= o in the other two is close 
to the ideal hydrodynamic limit at least until 1.5 fm/c. Later the expansion becomes three- 
dimensional and gluons in the outer regions cease to interact when the energy density de- 
creases under the critical value of lGeV/fm 3 . Therefore, the decrease in the transverse 
energy slows so that the final value of dET/dy\ y= Q is about 650 GeV, which is comparable 
with RHIC data 3- 

Although the interaction details do not matter for kinetic equilibration and pressure 
buildup, they do for chemical equilibration. Elastic collisions conserve the absolute par- 
ticle number and do not contribute to chemical equilibration, whereas multiplication and 
annihilation processes can drive systems toward chemical equilibrium. For the gluon evo- 
lution in central Au+Au collisions the initial free streaming (or the off-shell propagation) 
undersaturates the gluons (see Fig. [T5|) . For the pQCD bremsstrahlung processes chemical 
equilibrium is achieved by producing gluons. This leads to a larger Debye screening mass 
than that for elastic-only collisions (see Fig. [5]). Therefore, the elastic pQCD cross section 
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obtained from the simulation including the pQCD bremsstrahlung processes is the smallest 
(see Fig. Hj). 



B. Mean free path, mean transport path, and relaxation time 

Figure [10] shows the mean free path and the mean transport path of gluons. The mean 
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FIG. 10: Mean free path and mean transport path of gluons for both elastic and inelastic pQCD 
scatterings (solid curves), for elastic-only pQCD scatterings (dashed curves) and for elastic-only 
collisions with artificially large cross sections (dotted curves). 

transport path is defined as the inverse of the total transport collision rate and it is the path 
needed for gluons to reach isotropy in momentum space in a static medium. From Fig. [10] 
we see that the mean paths are all small early on and increase throughout the course of the 
expansion when the system becomes dilute. Comparing the mean free paths there is little 
difference for the processes with and without pQCD bremsstrahlung; however, the mean 
free path is much smaller when artificially large cross sections are considered. We also see 
that the mean transport path is larger than the mean free path for elastic-only collisions, 
because elastic pQCD collisions have small-angle scatterings and, therefore, do not isotropize 
the momentum efficiently. However, when pQCD bremsstrahlung processes are included the 
mean transport path and the mean free path are quite similar, so their kinetic equilibration 
is efficient. 

The relaxation rate of momentum isotropization l/9 re \{t) is calculated directly from Fig. 
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[6] using Eq. 



and is shown in Fig. [IT] in comparison with the total transport rate 



-^drift + -^22 + -^23 + ^32- According to Eq. ffl~5]) they should be identical, which is indeed 
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FIG. 11: Relaxation rate of momentum isotropization compared with the total transport rate 
for both elastic and inelastic pQCD scatterings (thin and thick solid curves) and for elastic-only 
pQCD scatterings (thin and thick dashed curves). 

seen within the numerical uncertainty. This indicates that the transport rates were correctly 
extracted. For the first 2 fm/c of the gluon evolution for pQCD bremsstrahlung processes 
the time scale of momentum isotropization is 1 — 2 fm/c, which is about a factor of 5 times 
larger than the mean free path (see Fig. [TOl . 



C. Collision rate, transport collision rate, and transport cross section 

Here we compare the collision rates, the transport collision rates and the estimates using 
the transport cross sections with each other concentrating on the results from the simulation 
with both elastic and inelastic collisions. 

Assuming Eq. (l28j) . the transport collision rates are directly proportional to the transport 
cross sections [see Eqs. ( H6l) . ( 1471) . and (HHl) ]. which can be directly linked to the collision an- 
gle distribution. To see how they differ from the true transport collision rates we compare the 
transport collision rates, n(v re \<j tr ) and the collision rates in Fig. [121 Multiplication factors 
according to Eqs. (l4~6l) and (148]) allow for more convenient comparisons. If the assumption 
( T28|) is realistic, the curves according to the assumption and those for the true transport 
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FIG. 12: Collision rate, transport collision rate, and n{v Te \a tv ) of the various collision processes. 
The results are obtained from the simulation, including both elastic and inelastic pQCD-scattering 
processes, and are depicted by the solid, dashed, and dotted curves, respectively. 

collision rate are identical. If additionally the collision-angle distribution is isotropic, all the 
curves in each case lie on top of each other. The fugacity X g in Eq. ( T4"6i) is calculated by 
X g = n/n cq , where n eq = 16 T 3 /tt 2 is the gluon density at thermal equilibrium at temper- 
ature T = e/3n. The gluon density n and energy density e are extracted from the parton 
cascade. 

We first examine the rates for the elastic-scattering processes shown on the left in Fig. 
[T2l where there is only a small difference between the true transport collision rate and the 
reduced rate related to the transport cross section. The difference comes from the Lorentz 
boost from the center-of-mass frame to the laboratory frame. Furthemore, we see that the 
transport collision rate is much smaller than the collision rate, which is again due to the 
fact that the pQCD gg — > gg scatterings are small-angle scatterings and are not efficient for 
momentum isotropization. 

In the middle the rates for the gg — > ggg bremsstrahlung processes are shown where 
little difference is seen. The transport collision rate divided by the kinematic factor 3/2 is 
the largest rate, especially over the collision rate. The kinematic factor for the assumption 
fl28|) is exactly 3/2. However, in general it is only approximately equal to 3/2, because the 
decomposition (135!) for R^ 3 [see W%1 in Eq. (1331) ] is an approximation. The real kinematic 
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factor defined as A may be larger than 3/2, which would lower R^/A below R23. Even 
though the difference between |i?23 an d R23 is small, which indicates that the collision- 
angle distribution in gg — > ggg collisions is nearly isotropic. The same is also seen on the 
right in Fig. [T2J where the rates for ggg — > gg are nearly identical. 

Figure [12] shows that the reduced transport collision rates related to the transport cross 
sections do not differ very much from the derived transport collision rates. Generally, the 
Lorentz boost from the individual center-of-mass frame to the laboratory frame does not lead 
to a big effect on momentum isotropization. This is nontrivial. However, it provides a basis 
to understand thermalization within multiparticle reactions. Additionally, the transport 
collision rates derivation helps to obtain the kinematic factors in Eq. (fl6i) . which are essential 
in quantitative analyses but typically ignored in the literature 

Bremsstrahlung processes gg <-> ggg are suppressed by the LPM effect, which occurs 
when a parton undergoes multiple scatters with radiated gluons through a QCD medium 
(originally photons in the QED medium). The interference of radiated gluons leads to 
suppression of radiation of gluons with modes (w,k), where w and k denote the gluon's 
energy and momentum. Heuristically, there is no suppression for gluons with a formation 
time t = wjk\ smaller than the mean free path. This is called the Bethe-Heitler limit, where 
the gluon radiation induced at different space-time points in the course of the propagation 
of a parton can be considered as independent events. Events within the Bethe-Heitler 
regime are included in BAMPS. Other gluon modes radiation with coherent suppression 
completely drops out, which is the reason for the G function in the matrix element in Eq. 
(J1J). Including these events speeds up thermalization; however, implementing the coherent 
effect into a transport model where the Boltzmann equation is solved remains a challenge. 

Without implementing the LPM effect as a strict low momentum cutoff, the matrix 
element for gg «-> ggg in Eq. fll]) is dominated by collinear bremsstrahlung, although it is 
suppressed by the Debye screening mass as an infrared cutoff. Therefore, the larger collision 
angle from gg <-> ggg processes in comparison to elastic scatterings originates from the 
present implementation of the LPM effect. 

Because the angle of the radiated gluon relative to the collision axis 6 is related to the 
momentum rapidity y by cos# = tanhy, which leads to coshy =1/ sin#, the effect of the 
Q function in the matrix element (j3J), Q(k^A g — coshy), on the angular distribution of the 
radiated gluon can be understood. For small transverse momentum k±, which the radiation 
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favors, the rapidity y is small due to the G function if the mean free path A g is small. This 
leads to large- angle radiation. The larger A g , the more small- angle bremsstrahlung (with 
large y) occurs. 

The function results in a cutoff in the radiated gluon phase space. The corresponding 
total cross section is found by integrating the matrix element 

(49) 

where g^. = <?i_/ s ; &± — A 9 — ^-gV^! = m l>/ s > is the angle between kj_ and q_i_, 

and if is a function of q±, k±, y, and <\>. H (found in Appendix D of Ref. jlj]) appears after 
the integral over the radiating gluon momentum. In Eq. (j4"9~l) the lower cutoff for k± and the 
upper (lower) cutoff ±y m for y can be seen, where y m is the minimum among arcosh(fcj_A 9 ) 
and arcosh(l/2fcjJ (see Appendix D of Ref. Q]). To obtain the radiated gluon angular 
distribution, which depends on m 2 D and A g , one has to integrate over k±, <j) and q± in Eq. 
(j39]). This is already done in Ref. jlj]. The radiated gluon angular distribution and the 
distributions of the other two gluons were depicted in Fig. 49 in Ref. [5], where m 2 D = 0.05 
and A g = 4. The distributions are nearly isotropic. From the present BAMPS calculation 
fh 2 D ~ 0.1 and A g « 3 when pQCD bremsstrahlung is included. The value of fh 2 D is almost 
identical with the equilibrium value m 2 D /(s) = Aa s /3iT = 0.13 for a s = 0.3. The smaller 
mjj found in Ref. [5] is due to the slower chemical equilibration, because the initial system 
(using p = 2 GeV) is more dilute than that used in this article (using p = 1.4 GeV). For 
larger fn 2 D and smaller A 9 large-angle scatterings for gg — > ggg are favored. 

Because the radiation is dominated by l/kj_ we simplify the matrix element (J4]) by elim- 
inating the collinear term l/[(kj_ — q_i_) 2 + ^d] ~ V^i) to see the effect of the LPM 
suppression (Bethe-Heitler regime) on the radiated gluon angular distribution. Then the 
radiation can be factorized 



1/2 _ rVm \ /-arcosh^/ A 9 /2 ;-l/2coshj/ _ J 

_ dk± / dyy- = / — / _ cf A;_l =- 

1/A 3 J —ym rC_\_ J— arcoshA/Ag/2 Jcoshy/A g 

-arcosh v /A 9 /2 2 COsh y J-y/l-2/K g 1 — U 2 

where u = cos 8. The integrand approximately represents the radiated gluon angular dis- 



tribution, which is bounded by ±y 1 — 2/A 9 . Figure [TBI shows the distribution for various 
A g = Ag^fs ~ A g T. The dist r ib ut ions are normalized and symmetric in cos9. The angular 
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FIG. 13: Radiated gluon angular distribution in the center-of-mass frame. 



distribution is peaked in the forward direction only for large A g as in elastic scatterings. In 
BAMPS we find A g « 3. Therefore, the radiated gluon angular distribution according to 
Eq. f|4"9]) is similar to A 9 t/s = 3 in Fig. [131 which indicates that large-angle radiation is 
favored. 



D. Dependence of the transport rate on the definition of Q 

We have already mentioned that the transport rates depend on the definition of the degree 
of momentum isotropy Q. In the previous subsections the numerical results for the transport 
rates with Q = (p 2 z /E 2 ) were shown. But the dependence of the numerical results of the 
transport rates on Q (specially when it is set to Q = (\p z \/E)) remains to be calculated. 

Reasonable definitions of Q must consider some kind of average of the momentum spec- 
tra shown in Fig. [2j thus, momentum isotropization time scales obtained from different 
prescriptions cannot differ much from each other. Because the inverse of the total transport 
rate is the momentum isotropization time scale, we do not expect any significant dependence 
of the transport rate on the definition of Q. 

Here we first compare the transport collision rates in the reduced formulas (l4*6l) for Q = 
(p 2 /E 2 ) with those in Eq. fj47j) for Q = (\p z \/E), where transport cross sections are defined 
by ([1]) and ([2]), respectively. We have already shown that if the collsion angle is isotropically 
distributed, the transport collision rates of a certain type of scattering processes are the 



34 



same, regardless of the definition of Q. For small-angle scatterings one has sin 2 9 ~ 9 2 
and 1 — cos9 « # 2 /2, and the transport collision rates for Q = (p 2 z /E 2 ) are a factor of 1.5 
larger than those when Q = (\p z \/E). For large-angle scatterings sin 2 f3 rs 1 — cos# and 
then the transport collision rates for Q = (p 2 z /E 2 ) are in turn a factor of 3/4 smaller than 
those when Q = (\p z \/E). The maximal relative difference amounts to 50%. Because pQCD 
bremsstrahlung is the dominant process in kinetic equilibration and the collision angle for 
that process is roughly isotropic due to the LPM cutoff, the difference in the transport 
collision rates due to different Qs should be minimal. 

Figure [TH shows the numerical results for the transport rates with the Qs defined above. 
The differences are only small, which means that the transport rates are not dependent on 
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FIG. 14: Transport rates calculated with Q = (\p z \/E) (solid curves) and Q = {p 2 z /E 2 } (dotted 
curves) . 

our choice in Q. 





E. Relaxation time r re i in the relaxation time approximation 

The collision term of the Boltzmann equation f)13p can be written as 

C( X ,p) = f ^ P \- f ^ P \ (51) 
r ie \{x,p) 

which describes the relaxation of the particle-density function f(x,p) by using a space-time 
and momentum-dependent relaxation time T re \(x,p), where r Te i(x,p) is a functional of f(x,p). 
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The underlying approximation in the so-called relaxation time ansatz is the ignorance of 
the momenum dependence of the relaxation time, i.e., r re i(x,p) ~ T ie i(x). Thus, r re i(x) 



gives the time sea! 



in 



e of the overall equilibration in absence of particle drift, which was used 



0, Q, Q Q, Q, , [ijj] to calculate the time scale of thermalization within various 
dynamical scenarios of the expansion. It is crucial to see whether r re i in the relaxation time 
approximation is equivalent to the mean transport path, because the latter determines the 
momentum isotropization time scale in a static system. 

For kinetic equilibration we insert (15T!) into the time derivative of the momentum 
isotropization ( [121) and obtain [by dividing (Q eq — Q)} 

1 f d 3 p p 2 z /eq-/ f d 3 p f cq -f\ 



R-22 + -^23 + -#32 - 7, _ / /<-, Q(t) I ., 

Q eq -Q\nJ (2ir) 6 E 2 r rc i nJ (2-Ky r rcl 

1 j ^eq Qcq ?l Q ^ ^eq ^ 



Q eq - Q \ n {T iel ) k n (r rol ) c y 

where (r re i)fc and (r re i) c are defined as averaged quantities over the momentum, and the 
index k denotes kinetic equilibration due to the convolution of angles (p 2 z /E 2 ) in the first 
integration, whereas c denotes chemical equilibration. 

(7rei)c can a ls° be calculated by integrating the collision term of the Boltzmann equation 
over the momentum 

' PP (C 22 + C 23 + C 32 ) = f LzzL = ^ , (53) 



which is a simple ansatz for the relaxation time 13| . The left-hand side of Eq. (153]) is equal 
to n(R 2 s/2 — -R32/3) if the explicit formulas of the collision terms and the definition of the 
collision rates are applied. We then obtain 

where the gluon fugacity is X g = n/n eq . 

Assuming that the relaxation time is independent of the momentum, (r re i)fc and (r re i) c 
become equal and one gets from Eq. f|52|) 

1 

\ E>tr 1 ptr 1 ptr 



(Trel)fe - . R tr , Dtr , Dtr ' ( 55 ) 
\ g n 22 i- -n-23 "r - n -32 



In chemical equilibrium (A fl = 1) the relaxation time is equal to the inverse of the total 
transport collision rate or the mean transport path. According to the relaxation time ap- 
proximation the right-hand sides of Eqs. (1541) and (1551) should be equal. However, it is not 
clear. 
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Without assuming the relaxation time ansatz we can also calculate the (r re i) c and (r re i)fe 
using Eqs. (154"]) and (1521) . because all the collision rates and the transport collision rates 
are known from numerical simulations. If the two "relaxation times" differ much, one can 
conclude that T Te \(x,p) is strongly momentum dependent and cannot serve as a global quan- 
tity to determine the overall gluon thermalization time scale in ultrarelativistic heavy-ion 
collisions. 

Before we calculate (r ro i) c and (r re i}ifc, we need the equilibrium particle density function 
f eq (x,p). Because we neglect quantum effects like gluon enhancement, f cq (x,p) = ve~ E l T at 
the center of the collision, where v = 16 is the degeneracy of gluons. The temperature T can 
be found using e eq = e, which stems from energy conservation in sudden thermalization 
However, the current particle density function could have an exponential shape, / = X g f eq , 
if the kinetic equilibration progressed quicker than chemical equilibration. In this case one 
obtains n = X g n eq as well as e = A 9 e eq and the temperature is then T = e eq /3n eq = e/3n. 
This temperature is larger (or smaller) than the previously defined temperature, if X g is 
smaller (or larger) than 1. The difference in these two local temperatures leads to the 
difference in e eq , n eq , and X g . Letting S denote sudden thermalization and E thermalization 
that follows an exponential behavior, it is easy to verify that 

*f = (^) 4 =^> (56) 

where they differ by a power of 4. The time evolution of the gluon fugacities A^ and 
X g obtained from BAMPS, including pQCD bremsstrahlung processes is shown in Fig. [T5l 
Although the system of minijets is initially slightly oversaturated, it becomes undersaturated 
due to a short period of (quasi-)free streaming. The reason is obvious from Eq. ( J56l) . when 
n as well as e decreases as 1/t in free streaming. The decrease of A^ is roughly a factor of 4 
stronger than that of A^ . Whereas Xg increases and relaxes to 1 later on, which indicates the 
ongoing chemical equilibration, A^ ~ 1 throughout the entire expansion, which implies that 
the system is in chemical equilibrium. The difference between A^ and A^ can be understood 
according to Eq. ([56]), so (A^ - 1) « 4(Aj - 1) for \X s g - 1| < 1. 

Physically, fugacity is a quantity that balances particle production and annihilation. 
Therefore, the ratio of the annihilation rate -R32 to the production rate -R23 can serve as a 
quantitative measure of fugacity, so, 2i? 32 /3i?23 is shown in Fig. [TS] 2R 32 /3R 2 3 agrees well 
with Xg , which implies that A^ is an appropriate choice for the fugacity in this example. 
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FIG. 15: Gluon fugacity. The solid, dashed, and dotted curve depict, respectively, \g, \g, and 
2i?32/3i?23 obtained from the simulation, including elastic and inelastic pQCD bremsstrahlung 
processes. 

The exact momentum averaged "relaxation times" (r re i)fc and (r re i) c according to Eqs. 
( 1521) and (1541) are shown in Fig. [16] by various gluon fugacities. The mean transport path 
V(-^22 + -^23 + -^32) i s a l so depicted for comparison. Except for (r re i)fc with A^ all "relaxation 
times" are considerably larger than the mean transport path. When comparing (r re i)fc to 
( r rei) c with the same fugacities, one finds large differences, especially for A^, where both 
"relaxation times" deviate by a factor of 4 to 6. This implies that the relaxation time 
r re \(x,p) in Eq. (l5Tj) indeed has a strong dependence on the momentum. Therefore, the 
applicability of the T re \(x,p) ~ T re i(x) in studying gluon thermalization in heavy-ion collisions 
is questionable. 



VII. SUMMARY 



Employing our recently developed parton cascade BAMPS and including inelastic pQCD 
bremsstrahlung processes we have introduced and calculated the transport rate of gluon 
drift and the transport collision rates of various scattering processes within relativistic ki- 
netic theory. We try to explain the observed fast equilibration of gluons within BAMPS in 
theoretical terms. 

We have shown that the derived transport rate of a certain process, -Rd rift , -^22 > -^23 > or -^32 > 
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FIG. 16: Relaxation time. The thick (thin) solid curve depicts the momentum averaged "relaxation 
time" (r re i)fc ((T rc i) c ) using the fugacity A^. The thick (thin) dashed curve depicts (r re i)fc ((T re i) c ) 
using the fugacity A^. The dotted curve shows again the mean transport path (see Fig. ITU]) . 
Results are obtained from the simulation, including both elastic and inelastic pQCD-scattering 
processes. 

determines exactly the contribution of the process to the denned momentum isotropization 
with Q = (p 2 z /E 2 ) (or Q = (\p z \/E)). The total transport collision rate, R\\ + R% + 
i?32, definitively describes momentum isotropization, whereas the change of gluon drift is 
a consequence of collision processes and -Rjf ri f t is negative in an expanding system. The 
inverse of the total transport rate, l/(i?^. ift + R\\ + + -R32), gives the exact time scale 
of momentum isotropization 9 ie \, and is about 1 fm/c from BAMPS for the gluon matter 
produced at RHIC. It is also shown that the calculated transport rates are independent on 
the definition of the degree of momentum isotropy Q. 

The inclusion of quarks into BAMPS is straightforward, but it is not yet completed. In 
the presence of quarks the Debye screening mass will be slightly larger, which leads to a 
decrease in the cross sections for gg — > gg and gg — > ggg scatterings. This slightly slows 
thermalization. However, further kinetic processes like q + g^q + g,q + g^q + g + g 
and q + q^q + q + g will speed up gluon thermalization. In addition, the effect may be 
small, because the initial quark density is tiny (30%) in comparison to the gluon density. 
Therefore, quark thermalization may be quite slow, but the study of it is still in progress. 
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We also derived the reduced transport collision rates related to the transport cross sec- 
tions. They are only exact for the special case when the center-of-mass frame of individual 
collisions coincides with the laboratory frame where the medium is observed. The devia- 
tions from the exact transport collision rates stem from the effects of Lorentz boosts from the 
center-of-mass frame to the laboratory frame. The numerical results show that the reduced 
rates differ little from the exact ones for the evolution of gluons in relativistic heavy-ion 
collisions. Lorentz boosts do not seem to have a major effect on momentum isotropization. 
Nevertheless, the derivation of the transport collision rate helps to obtain the correct kine- 
matical factors in the reduced rates summarized in Eqs. (j4"6]) and (T4T1) . which have been 
typically ignored in the literature. For instance, our analyses showed that a 2 — > iV (N > 2) 
production process is about a factor of (N — 2)/2 more efficient for momentum isotropization 
than its back-reaction or an elastic-scattering process. 

Using the numerical results of the transport collision rates for the various scattering pro- 
cesses we have investigated the importance of including the pQCD bremsstrahlung processes 
in thermalization. The inclusion of the pQCD bremsstrahlung processes and their back re- 
actions, as implemented in BAMPS, increases the efficiency for thermalization by a factor 
of 5. Overall kinetic equilibration and pressure buildup have a time scale of about 1 fm/c. 
The large efficiency stems partly from the increase in the particle number for the final state 
of gg — > ggg collisions but mainly from the almost isotropic angular distribution in the 
bremsstrahlung process due to the effective implementation of LPM suppression, which still 
needs to be further developed. 

Additionally, we have calculated the momentum averaged "relaxation times" with various 
gluon fugacities and they differ significantly from each other. This indicates a strong mo- 
mentum dependence of the gluon relaxation time r re i in heavy-ion collisions. Thus, using the 
standard relaxation time approximation of full kinetic Boltzmann processes is questionable. 
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APPENDIX A: DEGREE OF MOMENTUM ISOTROPY IN THE CENTRAL 
REGION 



The central region is described by a cylinder with a radius of r& = 1.5 fm and a longitudinal 
extension of 2Z b . The longitudinal boundary Z b = (tann?^) t with r] b = 0.2 increases linearly 
with time. Within the central region the degree of momentum isotropy is defined by 



Q(t) 



1 r d 3 p p 2 z 1 r'b r 2 ^ 
nlj^ffE^VJo rr Jo ' 

1 f d 3 p pi 2 f r b r2n 

drr 



where 



n(t) 



n J (2tt) 3 E 2 V Jo 
d 3 p 2 rn 



dzf(x,t,p) 
dzf(x,t,p) , 



, drr 
(2tt) 3 V Jo jm 



2tt 



dzf(x,t,p) 



(Al) 



(A2) 



V = 2irr 2 Zb is the volume of the central region. The second equation in Eq. ( lAlj) arises 
because of the symmetry of f(x,t,p) under the x —* — x exchange. In the limit — > and 
r/ b ^ one has 

\ rr b r2ir rZ b 

— drr d<p dz f(x,t,p) f(x,t,p)\ s=0 , (A3) 
V Jo Jo J-z b 

which is the definition of the degree of the local momentum isotropy in this limit [see Eq. 
( TTUj) in Sec. |Vj. The transport rates in this limit were already given in Sec. FVl , 
Taking the time derivative of Q(t) yields 



o 



dz — + taiahr] b f(x ± ,Z b ,t,p) 



z b Qf 

dz — + tanh rj b f(x±, Z b , t, p) 
o at 



■ (A4) 



The second term in the brackets comes from the time derivative of the boundary Z b and can 
be rewritten as 

tanhr] b f(x_ L ,Z b ,t,p)=[ dz ta ^ 17 ^ 6 f(x±, Z b , t,p) = / dz\f(x ± ,Z b ,t,p). (A5) 

Jo Z b Jo t 

The Taylor expansion of f(x±, Z b ,t,p) at x to the first order yields 



f(x ± , Z b , t,p) ~ f(x, t,p) + {Zb - z) 



(A6) 
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and we then obtain 
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(A7) 



(2tt) 3 V Jo 

The Oth-order contributions in Eq. (IA4I) cancel due to the definition of Q(t). (Z b — z)/t 
expresses the relative velocity of the boundary slice at Z b to the slice at z where particles 
are sitting. The second term in the brackets in Eq. (1A7f) appears due to the increasing 
longitudinal boundary of the central region, and it becomes smaller when Z b — >• (or 
r]b — > 0). According to the Boltzmann equation ffl3|) the expression in the brackets can 
be written as 



df (Z b - z) df 
dt t dz 



Vx_df_ _ Py_df_ 
E dx E dy 



E 



z\ df 



^ + C 22 + C 23 + C 32 . (A8) 



The term p z / E — (Z b — z)/t implies that only particles with longitudinal velocity p z j E larger 
than the relative velocity of the boundary slice (Z b — z)/t can drift out of the central region. 
This will be taken into account when calculating the transport rate of particle drift within 
the central region. The evaluations of the transport collision rates R 22 , R^, and -R32 are 
more straightforward. One only needs to replace the expressions derived in the limit r b — ► 
and T] b — > 0, which are already given in Sec. [V] by 



R 



tr 



1 

V Jo 



drr 



2- 



dzR) 



(A9) 



where i = 22, 23, or 32. 



[1] M. Gyulassy and L. McLerran, Nucl. Phys. A 750, 30 (2005). 

[2] P. Huovinen, P.F. Kolb, U.W. Heinz, P.V. Ruuskanen, and S.A. Voloshin, Phys. Lett. B 503, 
58 (2001). 

[3] S.S. Adler et al. (PHENIX Collaboration), Phys. Rev. Lett. 91, 182301 (2003). 

[4] J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 92, 052302 (2004). 

[5] Z. Xu and C. Greiner, Phys. Rev. C 71, 064901 (2005). 

[6] Z. Xu and C. Greiner, Nucl. Phys. A 774, 787 (2006). 

[7] P. Danielewicz and M. Gyulassy, Phys. Rev. D 31, 53 (1985). 

42 



[8] M.H. Thoma, Phys. Rev. D 49, 451 (1994). 

[9] G. Baym, Phys. Lett. B 138, 18 (1984). 

[10] A. Hosoya and K. Kajantie, Nucl. Phys. B 250, 666 (1985). 

[11] S. Gavin, Nucl. Phys. B 351, 561 (1991). 

[12] H. Heiselberg and Xin-Nian Wang, Phys. Rev. C 53, 1892 (1996). 

[13] S.M.H. Wong, Nucl. Phys. A 607, 442 (1996). 

[14] A. Dumitru and M. Gyulassy, Phys. Lett. B 494, 215 (2000). 

[15] J. Serreau and D. Schiff, J. High Energy Phys. 0111, 039 (2001). 

[16] P. Danielewicz and G.F. Bertsch, Nucl. Phys. A 533, 712 (1991). 

[17] A. Lang, H. Babovsky, W. Cassing, U. Mosel, H.G. Reusch, and K. Weber, J. Comp. Phys. 

106, 391 (1993). 

[18] W. Cassing, Nucl. Phys. A 700, 618 (2002). 

[19] D. Molnar, Nucl. Phys. A 661, 205 (1999). 

[20] J.F. Gunion and G. Bertsch, Phys. Rev. D 25, 746 (1982). 

[21] T.S. Biro, E. van Doom, B. Miiller, M.H. Thoma, and X.-N. Wang, Phys. Rev. C 48, 1275 
(1993). 

[22] T.S. Biro, B. Miiller, and X.-N. Wang, Phys. Lett. B 283, 171 (1992). 

[23] S. Mrowczyhski, Phys. Lett. B 314, 118 (1993); Phys. Rev. C 49, 2191 (1994); Phys. Lett. B 

393, 26 (1997); Acta. Phys. Pol. B 37, 427 (2006). 

[24] P. Arnold, J. Lenaghan, and G.D. Moore, J. High Energy Phys. 0308, 002 (2003). 

[25] P. Arnold, J. Lenaghan, G.D. Moore, and L.G. Yaffe, Phys. Rev. Lett. 94, 072302 (2005). 

[26] A. Rebhan, P. Romatschke, and M. Strickland, Phys. Rev. Lett. 94, 102303 (2005). 

[27] A. Dumitru and Y. Nara, Phys. Lett. B 621, 89 (2005). 

[28] B. Schenke, M. Strickland, C. Greiner, and M.H. Thoma, Phys. Rev. D 73, 125004 (2006). 

[29] K.J. Eskola, K. Kajantie, and J. Lindfors, Nucl. Phys. B 323, 37 (1989). 

[30] X.N. Wang and M. Gyulassy, Phys. Rev. D 44, 3501 (1991). 

[31] K.J. Eskola and X.N. Wang, Phys. Rev. D 49, 1284 (1994). 

[32] L.D. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 (1994); ibid. 49, 3352 (1994). 

[33] J. Adams et al. (STAR Collaboration), Phys. Rev. C 70, 054907 (2004). 

[34] J.D. Bjorken, Phys. Rev. D 27, 140 (1983). 

[35] H. Satz and D.K. Srivastava, Phys. Lett. B 475, 225 (2000). 

43 



[36] S.A. Bass, B. Miiller, and D.K. Srivastava, Phys. Lett. B 551, 277 (2003). 



44 



